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ABSTRACT 

We present a new Subaru/HiCIAO high-contrast 17-band polarized intensity (PI) image of a nearby 
transitional disk associated with TW Hydrae. The scattered light from the disk was detected from 0"2 
to T/5 (11 - 81 AU) and the PI image shows a clear axisymmetric depression in polarized intensity 
at ~ Of !4 (~ 20 AU) from the central star, similar to the ~ 80 AU gap previously reported from HST 
images. Azimuthal polarized intensity profile also shows the disk beyond 0'/2 is almost axisymmetric. 

We discuss two possible scenarios explaining the origin of the polarized intensity depression: 1) a 
gap structure may exist at ~ 20 AU from the central star because of shallow slope seen in the 
polarized intensity profile, and 2) grain growth may be occurring in the inner region of the disk. 
Multi-band observations at NIR and millimeter/sub-millimeter wavelengths play a complementary 
role in investigating dust opacity and may help reveal the origin of the gap more precisely. 

Keywords: planetary systems — protoplanetary disks — stars: pre-main sequence — stars: individual 
(TW Hydra) — techniques: polarimetric 


1. INTRODUCTION 

The structures of young circumstcllar disks can 
reveal key insights into the formation and evolu¬ 
tion of planetary systems. Complex structures such 
as gaps, holes, spiral arms, or asymmetric struc¬ 


tures have been imaged in a number of protoplane¬ 
tary disks at near-infrar ed (NIR) and mi l limeter/sub- 


millimeter wavelengths (Hashimoto et al.| 20111. 

2012; 

Andrews et al. 2011tlMuto et al.ll2012t Isella et al.l 

2013; 

Gradv et al.l 120131 iFukawaea et al.l 120131: iMavama et al. 
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I2012t iCurrie et al.l I2014H and they are predicted in 
theoretical simulations throughout the literatur e: 


planet embedde d in a viscous disk (IMuto et. al 


Zhu ct all l201li l. ph otoevaporation dClarke et al 


2010; 


2001 


Anu ct al.l I2U11II. p __ 

Owen ct al.i ; 2t)llaF10wenct al.ll2011b j. and grain growth 
dBirnstielet alJl2012hfl . 

The 7-10 Myr old (iDucourant et al.l I2014H pre-main 
sequence M2 dDebes et al.ll2013fl star TW Hydrae (here¬ 
after TW Hya) shows particularly strong evidence for 
a disk making a transition from a pre-planet build¬ 
ing mor phol ogy to a post-planet building structure (e.g. 
IStrom et al.l llOSl. TW Hya. w ith 0.55 ± 0.15 M 0 
in stellar mass dDebes et ahll2013fl an d more than 0.05 
Mq in disk mass dBergin et al.l 1201311 . is an ideal ex¬ 
ample of transitional di sks in that it is almost face-on 
(i = 7°; lOi et al.l 12004( 1 and beca use of its close dis - 
tance of 54 ± 6 pc (Hipparcos; Ivan Leeuwenl I2007H . 
is particularly accommodating to investigations of its 
geometrical structure with high spatial resolution and 
sensitivity. Previous studies at optical, near-IR, and 


ness extending to ~ 150 - 230 AU ( 

Krist et al.l 2000: 

Anal et al. 2004): Roberee et al.l 20051: 

Weinberger et al.1 

2002t Qi et al. 2004T 2006: Andrews e 

al.. 20TI7 High 


contrast direct imaging with the Hubble Space Telescope 
(HST) NICMOS coronagraph reveal a disk gap at 80 AU, 
indicating the evolution of dust composition or the pres¬ 
ence of an embedded protoplanetary system dDebes et al.l 
12013 1. The inner hole at ~ 4 AU is observationally pre¬ 
dicted from a flux deficit at ~ 2 - 20 pin in the spectral 
energy distribution (SED), suggesting a deficit of mate¬ 
rial left around the central star. Excess at ~ 20 - 60 pm 
in the SED is probably due to emission from the tru n- 
cated inner edge of the outer disk dCalvet et a 1.1121)0 2 i. 
The existence of the inner disk is also supported by in¬ 
terferometric observati ons and visibility model fitting of 
the 7 mm continuum dCalvet et al.l l2002t iHughes et alJ 
l2007fl and a recent modeling con cludes that the disk does 
not have a sharp edge at 4 AU dMenu et. alj|2014lf . 

In this Letter, we present a new 77-band polarized in¬ 
tensity (PI) image of TW Hya and the polarized inten¬ 
sity distribution ranging from 072 to 175 (corresponding 
to ~ 11 - 81 AU) of the disk, including the region closer 
to the central star than the region previously observed 
by HST. A polarized intensity depression at ~ 074 (~ 20 
AU), probably caused by a gap in the disk, was detected. 
Through investigating observed polarized intensity pro¬ 
file, we discuss the disk morphological structure and the 
origin of the depression. 


2. OBSERVATIONS AND DATA REDUCTION 

We carried out 77-band (1.6pm) linear polarimetric 
observations of TW Hya using the High Contrast In¬ 
strument for th e Subaru Next Gene ration Adaptive Op¬ 
tics (HiCIAO; ITamura et al.l I2006H with a dual-beam 
polarimeter at the Subaru 8.2 m telescope on March 
24th 2011 and January 3rd 2013. These observations 
were part of the Strategic Exp lorations of Ex oplanets 
and Disks with Subaru (SEEDS: ITamura 112009! ) project. 
Both of the observations were conducted in the same in¬ 
strument setting. The polarimetric observation mode ac¬ 
quires o-ray and e-ray simultaneously, and images a field 
of view 10" x 20"with a pixel scale of 9.50 mas/pixel. 

The exposures were performed at four position angles 


(P.A.s) of the half-wave plate, with a sequence of P.A. 
= 0°, 45°, 22.5°, and 67.5°, to measure the Stokes pa¬ 
rameters. The integration time per wave plate position 
was set by 10 - 30 s during observation to optimize in¬ 
strument according to on-site condition. Three coadd 
modes were used to save the t otal overhead time. The 
adaptive optics system fA0188: lHavano et al.ll200§l pro¬ 
vided mostly stable stellar point-spread function (PSF) 
with a FWHM of 0715 in the 77-band. The adaptive 
optics corrections in both observation dates functioned 
well and the qualities of resultant images of both obser¬ 
vations were almost the same. We strictly selected only 
good quality of data from among of them for producing 
the best image. Low quality images were removed prior 
to production of the final image, resulting in 13 good 
sets with a total integration time, for the PI image, of 
390 s. The data sets were taken with Angular Differen¬ 
tial Imaging (ADI) mode with a net field rotation of 12 
degrees. 

All unflagged data were corrected for field rota- 
tion and integrated. The standard imaging method 
(iHinklev et all I2009H was applied to the polarimet¬ 
ric data, including corrections for distortion and flat- 
held, using the Image Reduction and Analysis Facil¬ 
ity (IRA IL 1 ]) as i n pre v ious SEEDS publi cations (e.g. 
iHashimoto et al.l 120111 : IMuto ctahl I2012H . Correc¬ 
tion for instrumental polarization of the Nasmyth Hi¬ 
CIAO instrument was applied with the Mueller-matrices 
technique, that quantizes polarization characteristics 
Moos et all 120081 1. PI was calculated from \JQ 1 + U 2 , 
where Q and U are orthogonal linear polarizations, and 
the final polarized intensity profile was obtained after 
applying 10 x 10 pixels smoothing to the entire disk. 

3. RESULTS 

The PI image in 77-band of the transitional disk 
around TW Hya is presented in Figure Q] There is no 
discernible misalignment in the polarization vector be¬ 
yond r = 072 (which will be the subject of later studies 
in this series), and to be conservative, polarization mea¬ 
surements at radii more distant than 072 from the central 
star are regarded as more accurate. 

Shown in the PI image is a ring-like polarized intensity 
depression at r ~ 074. The observed polarized intensity 
profile shown in Figure[2ja), revealing a stair-like decline 
in polarized intensity, with increasing distance from the 
central star. We define three major step transitions as: 
Zone 1 (072 < r < 074), Zone 2 (074 < r < 078), and Zone 
3 (078 < r < 175), and the polarized intensity profile is 
fitted at each zone with a simple power law: oc r 1 . 
The fitting results in zones 1, 2, and 3 are shown in red, 
green, and blue, with slopes at each zone of -1.39, -0.33, 
and -2.65, respectively. The significant change in slope 
at ~ 074 is interpreted as a transition of some physical 
property in all azimuthal directions and is the subject of 
the remainder of this Letter. The change in slope ar ound 
078 has alr e ady b een discussed in lKrist et al.l (12000H and 
lAuai et ahl (| 2 00411 . 

Figure Ob) displays an azimuthal brightness profile at 

1 IRAF is distributed by National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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0"2, 0'/5, and 0"8 in radius. This figure shows the disk 
is almost axisymmetric, with constant polarized inten¬ 
sity in all azimuthal directions, although rms noise level 
increases at r = 0"2. 


T(r), has a form of T(r) oc r q , where q is a power index, 
c s oc r -<? / 2 , and ttk oc r ~ 3 / 2 , we obtain 

H s oc r^~ q ^ 2 . (6) 


4. DISCUSSION 

One of interesting results in the observation is that 
the polarized intensity profile has a stair-like shape. As 
shown in Figure [2ji, in Zone 1, the innermost part of the 
disk, the profile shows a shallow slope of approximately 
—1.4, where the slope is expected to be —2 in the case 
of a flared disk, or steeper in the case of a flat geometry 
(jWhitnev fc Hartmannlll992l ; llnoue et al.ll2008l ) . In Zone 
2, the profile appears almost flat, perhaps indicating that 
temperature increases with radius. In Zone 3, the polar¬ 
ized intensity profile has a steep value of approximately 
—2.65, which is darker than expected. The observed po¬ 
larized intensity profile reflects a radial distribution of 
small dust particles and dust property. In this section, 
we discuss three possible explanations for the radial po¬ 
larized intensity distribution observed in the disk of TW 
Hya. These include (a) an unusual radial distribution on 
the temperature, (b) the presence of a disk gap induced 
by the presence of a low-mass planet, and (c) a change 
in dust scale height due to grain growth. 


4.1. The Radial Temperature Gradient 

It is important to investigate the extent to which the 
shallow slopes observed in Zones 1 and 2 cou ld depend 
on temperature. Generally, llnoue et al. (200811 shows the 
observed intensity of the scattered light (J 0 bs) from a disk 
with no complex structures, such as gap or cavity, can 
be expressed by 

/ obs ~ (!) 

where L* is a stellar luminosity, r is a distance from a 
central star, and u(6) is a product of albedo and polar¬ 
ization (a phase function including polarization), 9 is a 
scattering angle, and /3 is grazing angle, and expressed 

= (2 ) 
dr r d r \ r J 


where H s is scattering surface. Assuming co[9) is con¬ 
stant in Equation m and considering J 0 b s oc r fc , we 
obtain 


d_ 

dr 



„i-fc 


(3) 


Thus, H s becomes 

H s oc r 3 ~ k . (4) 

According to customary practice and precedents in the 
past works, we assume that H s is linearly related to 
H g by some constant factor A in the case of smooth ly 
extended disk with no ga p dKcnvon fc Hartmannl 119871 
IChiang fc Goldreichl 119971) . and the relation can be de¬ 
scribed by 

H S ~A^ = A H s , (5) 

“fc 

where c s is the speed of sound, 12*, is Keplerian angu¬ 
lar velocity, and H g is a gas pressure scale height, re¬ 
spectively. Considering radial temperature distribution, 


For Zone 1, with radial polarized intensity dropping 
as r -1 ' 39 , Equation (JJ) implies H s proportional to r 161 
and q = —0.22. For Zone 2, with radial polarized inten¬ 
sity dropping as r -0 - 33 , H s is proportional to r 2 67 and 
q = —2.34. These indicate the temperature increases 
with increasing r, which is difficult to explain by con¬ 
sidering only stellar radiation or flared structure of the 
disk. One of the reasons for the unexpected result is that 
the gas scale height is not simply related to temperature. 
Such a result might be able to be explained when a gap 
exists because the gap generates a sudden change of sur¬ 
face density. Dynamical phenomena such as disk wind 
or turbulence could be cited as another possible cause, 
because the disk surface can be lifted by them, result¬ 
ing in a flat polarized intensity profile and no longer in 
hydrostatic equilibrium. If this is the case, wind mech¬ 
anism such as photoevaopration or magneto-rotational 
instability that suppress the inner region of the disk due 
to dead zone might occur. 

4.2. Planet-Induced Disk Surface Structure 

In this section, we consider the possibility that the 
changing slope of the observed polarized intensity dis¬ 
tribution may signal the presence of a planet within the 
disk. The slopes at each zone ca n be explained by a disk 
surfac e geom etry described by iJang-Condell fc Turned 
(|2012t 1201311 . In their works, the disk illumination 
and shadowing are self-consistently calculated consider¬ 
ing several disk orientations, using a static a-disk radia¬ 
tive tra nsfer model when a gap exists in the disk. Fig¬ 
ure 2 in lJang-Condell fc Turned (|2012l ) shows the surface 
brightness becomes brighter at the gap wall directly ir¬ 
radiated by a central star, and darker by shadowing and 
cooling behind the outer gap wall at the edge of gap. If 
this is the case, Zones 1, 2 and 3 probably correspond to a 
gap, an outer gap wall, and a shadowed region behind the 
outer gap wall, respectively. They also made radiative 
transfer calculations at several wavelengths and theoret¬ 
ically derived a radial surface brightness profile of a disk 
when a gap is formed by a 70 and 200 earth mass planet. 
The radial surface bright ness profile at a waveleng th o_f 1 
prn shown in Figure 6 in IJang-Condell fc Turned (120121 1 
increases just beyond the gap created by a planet, and its 
slope depends on planet mass, flatter in the case of small 
mass planet. Therefore, if the origin of the polarized 
intensity depression attributes to a planet, planet mass 
seems to be small since our observed polarized intensity 
profile shows a flat slope at the gap region. 

It is beyond the scope of this paper to investigate the 
detailed effects of planet parameters on the structure 
of the inner disk of TW Hya, but an extensive search 
of parameter-space (planet mass and orbital characteris¬ 
tics) in the future might lead to a suitable model for the 
observed brightness distribution. 

4.3. Possibility of Grain Growth 

The pictu re of the disk structure presented in Section 
14.11 and 14.21 does not necessarily mean that the gas disk 
has a gap at r ~ 20 AU. Dust grains are also believed 
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to settle to the mid-plane as their size increases. When 
the disk gas has a Gaussian distribution in the verti¬ 
cal direction with H g , dust grains also have a Gaussian 
distribution in the vertical direction, with a dust scale 
height H d that is determined by the balance between 
the diff usion by the gas and the coupling efficiency to 
the gas dDubrulle et al JUifObl) . Using the Stokes number 
St, which represents the coupling efficiency of grains to 
gas and is defined as the stopping time normalized by 
the orbital timescale, and the diffusion coefficient qd, 
which represents the diffusion efficiency of the gas in the 
vertical direction, the dust scale height can be written as 


//,] = H g max 1, 



(7) 


when St <C 1 (see iBrauer et all 120081 and references 
therein). As a first rough step to quickly quest to see 
a possibility whether grain growth occurs or not, we in¬ 
troduce the parameter s where the radial grain size dis¬ 
tribution has a form of a(r) oc r~ s . By using p in the 
form o f E g oc r~ p and applying turbulent viscosity a 
model (IShakura fe Sunvaevlll973l) and the Epstein law, 
then the Stokes number has a simple form of 


St ~ oc r p ~ s , (8) 

where p d is the density of a dust particle and assumed 
to be a constant throughout the disk. Substituting this 
into Equation 0 , we obtain 

H d oc iJ g r (s -P )/2 . (9) 

Since H g oc A 3 ~ q ^ 2 : we finally obtain the radial depen¬ 
dence of the dust scale height as 

H d oc r {3+s -P- q ' >/2 . (10) 


Our observation suggests that the dust scale height has 
a radial dependence of H d oc r 1 ' 61 in Zone 1 and H d oc 
r 2 ' 6 ' in Zone 2, by utilizing Equation 0 and assumption 
that H s is linearly related to H g described by Equation 
©• Note that H d and H s are not identical, however H s 
must be about equal to H d for the small dust grains, 
which are mainly traced by iJ-band scattering light in 
our observation. From Equation (fl(71) . we obtain s = 1.72 
in Zone 1 and s = 3.84 in Zone 2 in the dis k around TW 
Hya assuming that p = 1 and q = 0.5 dHuehes et alJ 
120081 ). This suggests that the dust has a larger size in 
the inner part of the disk, and the size varies by at least a 
factor of 10 between 10 AU (~ 0"2) and 40 AU (~ 0"8), 
where the depression and its vicinity are covered. 

It should be noted that since smaller grains have a 
larger scattering efficiency, the scattering plane has to 
represent the dust scale height of the smallest grains at 
the radial position. The PI image suggests the possible 
existence of a small-grain depleted region at r < 0"4, 
with the depleted range of width more than 10 AU wide. 
This may suggest that the larger grains are present at 
inner radii, but there is no sharp transition of the grain 
size distribution. We note that gas and large grains can 
be present in the inner region of the disk. Indeed, mm- 
/submm-observations give no evidence to date of the de¬ 
pletion of gas or dust continuum emission. Furthermore, 
the dust sedimentation near the disk mid-plane results in 


flat polarization intensity profile, that is it further pro¬ 
motes a flared structure in H d distribution. 

Estimates of the amount of 10 pm - 1 mm sized dust 
in the disk mid-plane, and up to ~ 60 AU, are limited 
by SMA and VLA data (see Fig. 9. iMenu et. all 12014 ). 
Our NIR observations trace dust at the disk surface be¬ 
tween 10 - 80 AU, and reflect mainly smaller grains 0.1 
- 10 pm in size. The observations are complementary to 
each other, revealing radially and vertically distributed 
dust at different size scales. Note that the absorption 
efficiency also depends on grain size, and thus influ¬ 
ences on the vertical profile via temperature distribution. 
However, the absorption efficiency does not significantly 
change with dust grain sizes between 0.1 10 pm in the 

case of ff-band wavelength dMivake fe Nakagawalll993lh 

The variation of fundamental dust properties may also 
affect the observed polarized intensity distribution. For 
instance, ui(6 ) in Equation ([1} may not be constant. 
iTakami et al.l (1201311 explored the scattering properties 
for different grain size, and calculate polarization degree 
and polarization intensity normalized by the incident 
stellar light as a function of scattering angle. They show 
the polarization degree decreases with increasing grain 
size when 2-7ra < A, where a is a radius of a dust grain 
and A is wavelength of emission, respectivel y, resultin g in 
a dep olarization of the emission (see also [Fischer et al.l 
11994T) . On the other hand, when 2™ > A, polarization 
and polarized intensity decreases less significantly even 
though a becomes larger. Since the dust grain size di¬ 
rectly affects the polarized intensity, a variable ui(8) is 
consistent with the effect of grain growth. 

The minimum grain size depends strongly on the frag¬ 
mentation efficiency of dust-dust collisions in the disk. If 
fragmentation of dust is efficient, the scattering surface 
of the disk should be dominated by the small fragments. 
However, our observational analysis suggests that frag¬ 
mentation is not so efficient, particularly in the inner 
part of the disk. One possible reason for the inefficiency 
of fragmentation is the icy aggregates have a larger crit¬ 
ical velocit y before complete disrup tion, than do rocky 
aggregates IWada et al.l (12009 ,. 20131 1. although the crit¬ 
ical velocity of fragmentation is still under debate. In 
our study, it is reasonable to assume that the grains are 
covered by ice, since the locations we have observed are 
expected to be outside the water snow line. 

Another possibility is that fragmentation is effective, 
but the tiny f ragmen ts are swept up by larger grains. 
[Birnstielet al. (2012b) investigated the size distribution 
as a result of the balance between fragmentation and 
sweep-up. To sweep up tiny grains, they suggest the 
diffusion coefficient a should be as low as «d = 10 -4 . In 
our case, the modeling therefore suggests the diffusion is 
generally inefficient throughout the disk. 

In tandem with NIR high-resolution observations, 
high-resolution and multi-band observations at mm- 
wavelength will aid an investigation into the grain opac¬ 
ity index /3. Such data will help us to clarify whether 
the grain growth is occurring at inner radii and under¬ 
stand the effects of any planet formation on dust char¬ 
acter within the gap. 
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Figure 1 . H -band PI image of TW Hya in the ADI mode. North 
is up and east is to the left. The dark filled circle at the center 
indicates a software mask with r = C//2. (A color version of this 
figure is available in the online journal.) 
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Figure 2. Polarized intensity profiles of the TW Hya disk at H 
band (left), and azimuthal brightness profile images at r = 0'.'2 
11 AU), 0/5 (~ 27 AU), and O'/8 (~ 43 AU), respectively (right). 
The error bars shows 1 cr. The numbers 1—3 in the left panel denote 
the three radial zones for which a separate polarized intensity slope 
is fit. 










































